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The developments in printing technologies allow fabrication of micron-size nano-layered delivery systems 
to personal specifications. In this study we fabricated layered polymer structures for drug-delivery into a 
microfluidic channel and aimed to interferometrically assure their topography and adherence to each other. 
We present a scanning white light interferometer (SWLI) method for quantitative assurance of the 
topography of the embedded structure. We determined rapidly in non- destructive manner the thickness and 
roughness of the structures and whether the printed layers containing polymers or/and active 
pharmaceutical ingredients (API) adhere to each other. This is crucial in order to have predetermined drug 
release profiles. We also demonstrate non-invasive measurement of a polymer structure in a microfluidic 
channel. It shown that traceable interferometric 3D microscopy is a viable technique for detailed structural 
quality assurance of layered drug-delivery systems. The approach can have impact and find use in a much 
broader setting within and outside life sciences. 

Technology allows developing and fabricating pharmaceuticals in a tailored and personalized manner to meet 
the needs of individuals and patient groups 1 " 3 . The challenge is to deliver the benefits of this work to patients. 
Personalized medicine relies on diagnostic/genetic testing and targeted drug delivery. In addition our vision 
of future personalized medicine also entails flexible fabrication to allow accurate dosing of medicines. 
Personalized medicines can be e.g. produced by a scalable technology based on precise and robust 3D printing. 
Real-time quality control of the functionality of these drug delivery systems would be important and desirable. 

A small, multi-layered, polymer structure can allow precisely controlled drug release in time and space with a 
carefully engineered release profile 4,5 . It can allow several active pharmaceutical ingredients to be released in 
temporal sequence. To assure the functionality of such a structure the thickness of the different layers should be 
tightly controlled, the layers need to adhere to each other, and the surface of the outermost layers must be intact 6 . 
The absolute measurement precision must be high because the layers are thin and because the drug release 
depends on layer thickness 7 and layer roughness 8 . The quality assurance method should be sterile and non- 
contacting to avoid contamination and altering the surface. 

Our proposed contribution to personalized medicine is a method that combines 3D printing, 3D scanning 
white light interferometry (SWLI), and microfluidics to enable engineering and control of drug release profiles. 
Each of these three primary methods already alone provides impact. 3D printing gives reproducible otherwise 
unobtainable release profiles (it also permits high spot density, small volume spots and microarrays) 9 , SWLI 
provides non-contact, real time and non-invasive structural inference that can be metrologically traceable. 
Microfluidics can provide a clinically relevant and controllable, small volume environment where temperature 
and (bio) chemistry can be rapidly changed 10 . 

In drug design and development one needs to control the properties of drug substances and additives in 
addition to enhancing and controlling their delivery and release 3 ' 6 . Printing drugs allow fabricating structured 
and tailored drug-delivery systems. One can use topology (e.g. thickness and surface roughness of the printed 
layers) and chemical concentration to engineer their release profile 6 . These systems can be produced to tight 
specifications. Since material can be deposited accurately (space and amount) printed drugs with predetermined 
release profiles can potentially carry low drug concentrations. Quality assurance of printed products remains a 
challenge 2 ' 4 . 

2D and 3D printing permits fabricating functional structures of metals, polymers, and biomaterials 1112 . 3D 
structures can be printed on a variety of surfaces with characteristic permeability, porosity, hydrophobicity/ 
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Figure 1 | Scanning White Light Interferometer set-up that compensates for the microfluidic channel roof (Left). The geometric thickness (d g ) of 
the polymer layer is A1-B2 since the refractive index of air = 1. The effective refractive index of the adhered layer is (B3-B2)/(A1-B2). The same approach 
can be used to determine the height of disbonds. The method provides 10-20 nm traceable resolution along the z-dimension. It can measure stacked 
layers that are static or moving 10 . It can determine a Sq parameter (ISO 25178) that potentially can be linked to Rq even inside a microfluidic channel 
without direct tactile access. 



hydrophilicity and surface energy. This allows controlling the prop- 
erties of the printed substances. Besides accurate patterning, printing 
on tailored functionalized substrates and multi-layer printing makes 
automated high-speed manufacture of complex structures possible. 
It opens up new avenues for tailoring physicochemical properties of 
organic substances. 

There are few non-contacting and non- destructive ways to mea- 
sure, in a traceable manner, surface roughness, thickness and adhe- 
sion of layers. Optical confocal microscopy provides sub -micron 
resolution along the z-direction (into the sample) but is disadvant- 
ageous for large area scanning. Moreover, it employs intense laser 
light, which may alter the mechanical and chemical properties of 
drug-laden polymers. Measuring inside a microfluidic channel is 
hard 13 . LED -based optical coherence tomography (OCT) even 
though a point-by-point method, fulfills many of the above require- 
ments but its precision in determining the thickness of thin layers 
may be insufficient 14 . Ultra-sound microscopy especially used in 
second harmonic mode 15 provides qualitative images but seldom 
allows one to quantitatively determine the strength of adhesion or 
the size of disbonds along the z-direction. Moreover it can't deal with 
gas-filled voids and usually requires acoustic couplants that may 
modify the sample. Terahertz imaging can infer about layered struc- 
tures but does not provide sufficient quantitative spatial precision 
along the z-direction 8 . Nano-scale x-ray tomography potentially 
offers sufficient precision but is impractical since it is slow and uses 
ionizing radiation 16 . None of these methods fulfill the requirements 
presented by the current case (personalized drugs produced at point - 
of-care). Non-contacting, rapid, and non-destructive SWLI allow the 
possibility to determine surface roughness, the geometric thickness 
of layered structures, the adherence between the layers, and the 
thickness of each individual layer 17 . Here we report the first imple- 
mentation of scanning white light interferometry to study these fea- 
tures of layered drug-delivery systems. 

Results 

Our custom-built SWLI optically compensates for the thickness of the 
roof in the microfluidic channel (Fig. 1) and measures the effective 
refractive index of the drug-laden structure (Fig. 2). The SWLI method 
provides 10-20 nm traceable resolution along the z-dimension. This 



potentially sterile technique requires no sample preparation. Using 
image stitching it allows, when necessary, high resolution also across 
a large surface (up to 200 X) 18 Moreover, it can determine surface 
roughness Sq (ISO 25178) for an embedded surface without direct 
tactile access. It can measure stacked layers that are static or 
moving 19 . 

SWLI measures the optical path difference between the reference 
and sample arm of the interferometer. When the sample is in air and 
we measure the step height - the single layer thickness - there are two 
interferences one for the base surface (Point Al, Fig. 1) and one for 
the top surface of the layer (Point B2, Fig. 1). The distance between 
these two interferences is equal to the layer thickness directly. When 
we measure through a well bonded layer there are also two interfer- 
ences, but the position of the second one is shifted (Point B3, Fig. 1), 
compared to the interference at the base surface (Point Al, Fig. 1) 
because the optical thickness is equal to the geometrical thickness 
multiplied by the refractive index of the layer. The shift of the second 
interference depends of the refractive index of the layer (case B in the 
figure). The ratio between geometrical and optical thicknesses is 
equal to the refractive index of the layer. In the case of an unbonded 
or missing part of the layer a third interference appears between the 
first and second one (case C in the figure) and the shift of the last 
interferogram (Point C6, Fig. 1) depends of the thickness of the 
nonbonded layer. 

To demonstrate the capacity of our proposed method a multi-layer 
drug-delivery system containing the vitamin B2, riboflavine sodium 
phosphate (RSP), is shown in Fig. 3a. Samples featuring these 
crossing lines of HPMC and RSP were placed on glass microscope 
cover slips. The image is 1 mm 2 . We prepared relevant test samples 
using a technique that simulated printing as described in 5 . The trans- 
parent HPMC lines were on average 3.98 |im thick and 5-10 mm 
wide whereas the RSP lines were on average 26 [im thick. The rough- 
ness of the HPMC (Sql) was 0.58 urn and RSP (Sq2) 0.46 urn 2 . 
Figure 3a shows a disbonded HPMC-RSP multilayer system featur- 
ing air filled voids 0.15-30 |im tall. Figure 3b illustrates both adhered 
and disbonded multi-layer systems of hydroxypropylcellulose 
(HPC) film on a PET (polyethylene terephthalate) surface. These 
images show high resolution top surface topography with roughness 
of 0.31 |im (Sql) and an apparently rougher second surface (light 
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Figure 2 | HPC-Theophylline laden layer on top of a HMPC layer residing on PET. Even without knowing the refractive index of the HPC-Theophylline 
layer the thickness can be determined (5 urn on the left 10 (am on the right). Sq on the indicated area is 0.157 (am. (Left) Measured interferences. (Inserts 
top right) Schematic of the samples. 
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Figure 3 | (a). Multilayer drug delivery system in the form of a crossing strip of printed riboflavine sodium phosphate (RSP) layer between two 
layers of hydroxypropyl methylcellulose (HPMC) film reconstructed from a SWLI measurement, (b). Structure of hydroxypropyl -cellulose (HPC) film 
on a PET surface revealing up to 13 (am tall delamination, 0.31 (am (Sql) and 0.20 (am (Sq2) surface roughness. (Left) Measured interferences, (Inserts 
top right) Schematic of the samples.We can determine surface topology with high resolution and we can determine presence of disbonds in multilayer 
structures. 
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Figure 4 | Left: (a) Interferograms and reconstructed SWLI image of printed ethylcellulose polymer film drug-delivery system imaged inside a 
microfluidic channel during 8 ul/min water discharge. Indicated: DDS surfaces (top surface and bottom surface). The single camera pixel data - 
interferograms, were obtained for the corresponding points on the layers. Using the method outlined in Figure 1 we determined the refractive index 
(n = 1.47 ± 0.02), the average thickness (44.0 (am) of the structure, and the average height of the disbond (5.7 (am). The Sq in the indicated area is 
0.30 (am. Right: (b) The imaged structure. 



dispersion by the top surface makes it appear rougher than it is). 
Single camera pixel data (Fig. 3a and 3b, leftmost) and a 3D surface 
plot of the polymer layer are shown (Fig. 3a and 3b, middle). The 
images were reconstructed pixel by pixel from the interferogram. The 
depth coordinate of each layer was extracted using Larkin's algo- 
rithm 20 . Fig. 2 shows how the geometric thickness of a structure 
can be determined even though its refractive index is not known in 
advance. The SWLI measures optical thickness d opt . The geometric 
thickness (d g ) of the polymer layer is A1-B3 since the refractive index 
of air = 1. The effective refractive index of the adhered layer is (B3 - 
B4)/(A1 - B3). The same approach can be used to determine both 
layer thickness and disbond height. Except for the top layer the depth 
coordinate was corrected using the refractive index. This index was 
obtained by relying on a discontinuity in the structure, see Fig. 3a. 
The interferograms from a point on the top layer and from a corres- 
ponding point on the bottom surface on that layer are separated on 
average by d opt = 6.23 |im whereas the discontinuity on bottom 
surface is known to be 2.25 um tall which indicates an estimated 
layer thickness of 3.98 urn. The height of the disbonds was deter- 
mined by assuming that the refractive index of the disbond was 1.00. 
Figure 4 shows a polymer structure measured inside a microfluidic 
channel during 8 ul/min (1.8 mm/s) water flow. It took a few min- 
utes to obtain the images in Fig. 3 and 4. These images were auto- 
matically analyzed by commercial and custom made software to 
construct and analyze 3D images of the cross -over point of the two 
crossing strips. The standard uncertainty for the z-coordinate in the 
construct was 0.045 urn. We assumed the refractive index of water to 
be 1.33. 

Discussion 

The proposed method rapidly determined the structure and presence 
of disbonds in the printed drug- delivery system. There are limita- 
tions to the method. A layer separation, of less than 0.5 um or more 
than a few millimeters is hard to image. When the number of inter- 
faces exceeds four it is hard to automatically extract the layers. The 
employed technique to derive the effective refractive index may not 
always work well. One example is high throughput analysis of sys- 
tems with different drugs, another example is if the drug concentra- 
tion exhibits a gradient inside a layer. Finally, the Sq value has not yet 
been linked in a traceable manner to a common surface roughness 
index such as Rq. 

We presented for the first time a practical method for traceable 
structural quality assurance of printed drug delivery systems. This 



result opens up a plethora of opportunities to induce quality assured 
fabrication in the field of tissue engineering, microfluidics, bio- 
MEMS, metamaterials, as well as nano- and micro-sized layered 
structures in addition to printed pharmaceuticals. 

This traceable quantitative method can provide a competitive edge 
to both research and production carried out in academia and the 
industry. In practice one can carry out high throughput label-free 
testing, one can optimize products structurally, and one can assess 
systems in wet and dry environments. The method has high spatial 
and temporal precision, and can be applied in a controlled micro- 
fluidic environment. It can allow structural and functional quality 
assurance for a broad range of applications in life sciences. 

Methods 

Scanning white light interferometry (SWLI). The samples were imaged using a 
custom-made SWLI instrument, Fig. 2. Briefly, the instrument uses a halogen lamp 
(Osram G4, driven at 6 Volts, 10 W), a standard 5X Michelson Nikon objective 
(Nikon CF IC Epi Plan TI; Japan) or 10 X Mirau objective (Nikon CF IC Epi Plan DI; 
Japan), and a piezo translator with 100 um travel (Physik Instrumente - type P- 
721. CDQ). The system magnification was 3.15 with the Michelson objective and 6.3 
with the Mirau objective. The instrument can determine profiles with ±15 nm 
accuracy along the z-direction. A detailed description of the VIS interferometer is 
given in 17 . In the receiving part we employed a 3.63 um pixel size black and white 
camera (Hamamatsu CI 1440 Orca Flash2.8, Hamamatsu City, Japan). 

We validated the SWLI for top surface imaging by inspecting a standard sample 
VLSI 1853 ± 2.3 A (model #SHS-1800 QC, VLSI Standards, Inc.). Traceability to the 
national SI standard is ensured using a transfer artifact. The influence of the roof of 
the channel was removed with a compensating plate identical to the 200 um thick 
cover glass (CLS2960246, Sigma-Aldrich, Corning) of the microfluidic channel 
inserted into the reference arm of the interferometer. 

Microfluidic channels. The cover glass and 1000 um thick microscope base glass 
were used (Menzel Microscope slides, ISO 8037/1) to construct the microfluidic 
channels. Inlet and outlet holes (2 mm in diameter) were drilled in the base glass. 
Plastic tubing (TubPEEK Blu 1/16 Z226661) was attached by UV cured glue 
(Thorlabs NOA81). The tubes were leveled to the base glass. A 2 mm wide channel 
was created using two-sided 50 um thick double-sided tape (Tesa) between the cover 
and base glass. The sample was placed in the middle of the channel. The cover glass 
was then attached. Ion-exchanged water was flowed through the channel with a pump 
(Alladdin-1000, World precision Instruments, Aston, Stevenage, UK). 

Sample preparation. There were four kinds of samples. The films were created by 
solvent casting using a manual pipette (BioHit ProLine + Mech, Finland) simulating 
printing of drug substances and polymers as described in 2 . All samples were dried in 
ambient conditions (21°C and 40% RH). The aqueous solutions for Figure 3 
membranes consisted either of 2 mg/ml riboflavin sodium phosphate (RSP, 
riboflavin 5' -monophosphate sodium salt, Ph. Eur., Fluka Analytical, Sigma-Aldrich, 
France) in 1% (w/w) hydroxypropyl methylcellulose (HPMC, Metolose 90SH-4000, 
Shin-Etsu, Tokyo, Japan) or pure 1% HPMC solution. The sample presented in 
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Figure 3a was a three-layer construct (HPMC/RSP-HPMC/HPMC, from top to 
bottom) prepared on a water impermeable polyethylene terephthalate (PET) film 
(Mylar® A, Dupont Teijin Films Europe, Luxembourg). The sample in Figure 3b was 
a two-layer construct with hydroxypropylcellulose (HPC, Klucel LF, Shin-Etsu, 
Tokyo, Japan) residing on a PET surface (HPC/PET, from top to bottom). The HPC 
film was prepared by deposition of 1% HPC solution on top of the PET surface. The 
solution to fabricate the film for the sample in Fig. 4 contained water insoluble 
ethylcellulose (EC, E8003, Sigma-Aldrich, USA; ethoxy content, 48-49.5%; viscosity, 
45 cP) and ethanol (>96.1%, Etax A, Altia OYj, Finland). The EC sample was placed 
in a microfluidic channel for measurement. Fig. 2 depicts a three-layer construct 
(theophylline-HPC/HPMC/PET, from top to bottom). The solutions to create these 
double-layer membranes contained 7 mg/ml theophylline in 3.3% (w/w) 
hydroxypropylcellulose (HPC, in ethanol) (top layer) and 1% HPMC (in water) 
(second layer). These two layers were made on top of a PET surface. Incisions through 
the double-layer membrane were made to remove a section of the system to allow 
determining the film thickness and its attachment to the PET-base in Fig. 2. 

Image reconstruction and processing. We scanned the samples with our custom- 
made SWLI instrument. Images were acquired with 68.75 nm vertical spacing. At 
each height we averaged 10 camera frames after which we recorded both the average 
image and the height coordinate output by the calibrated piezo capacitive sensor. 

Surfaces were reconstructed using Larkin's algorithm 7 on a higpass-filtered SWLI 
signal. This modification permitted us to do multi- surface feature extraction. Shortly: 
the single pixel camera intensity data vector (stacked pixels) was high-pass filtered 
with a second order derivative filter. The algorithm then identified and stored the 
maximum of the envelope of the interferences in the filtered vector. The interactive 
multi- interface extraction process first asked for the maximum number of interfaces 
to extract, and for the minimum permissible threshold value for envelope maxima 
and interface separation. For each 3D image the threshold value was found by an 
iterative procedure. (In Fig. 2 the HPC Theophyllin - HPMC interface was found by 
pooling 25 pixels prior to extracting envelope maxima. Extracted interfaces were 
filtered using first a 3 X 3 pixel filter and then a 7 X 7 pixel median filter along the x-y 
direction. Surface height data was corrected from optical height to geometric height 
using the method illustrated in Fig. 2. In practice, planes were fitted in least squares 
sense to the reference surfaces #1 and #3 in the above figure. The distance between 
these planes was determined analytically. Local surface fitting (second or third order 
in x and y) was used to patch small areas lacking data. These areas were identified and 
the size of the patch was decided manually post hoc. The final 3D image was tilt- 
corrected by requiring that the bottom most interface was horizontal. The Sq para- 
meter was determined from unfiltered surface data as outlined in IS025178. 
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